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The dynamics of observables which are matrices depending on h and taking values in classical phase 
space is defined retaining the terms up to the first order in h of the Moyal bracket. Within this semi- 
classical approach a first order lagrangian involving gauge fields is studied as a constrained hamiltonian 
system. This provides a systematic study of spin dynamics in the presence of non-abelian Berry gauge 
fields. We applied the method to various types of dynamical spin systems and clarified some persisting 
discussions. In particular employing the Berry gauge field which generates the Thomas precession, we 
calculated the force exerted on an electron in the external electric and magnetic fields. Moreover, a 
simple semiclassical formulation of the spin Hall effect is accomplished. 
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1 Introduction 

In [1] the intrinsic spin Hall effect was studied considering abelian and non-abelian Berry gauge 
fields [2] arising from the adiabatic transport. After this seminal work there has been a great ef- 
fort to employ Berry gauge fields to acquire a better understanding of spin dependent dynamics 
semiclassically |3]-|T7]. Although similar phenomena were treated letting coordinates and/or 
momenta be noncommuting, they appear to be disconnected. We would like to present a formu- 
lation which embraces these approaches. In our formulation keeping track of the semiclassical 
approximation is easy and interactions between different gauge fields can be introduced in a 
simple manner. 

To present our approach we need to recall the Weyl-Wigner-Groenewold-Moyal (WWGM) 
method of quantization as well as the Dirac formulation of constrained hamiltonian systems [19 

Quantum dynamics of particles without spin is usually provided by operators depending on 
the quantum phase space variables (p^, x^) satisfying the Heisenberg algebra: [p^, x^] = —ihS^^, 
[Pfj.jPu] = 0, = 0. However, there is an alternative approach due to WWGM where one 
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introduces symbols of operators and their star product: Observables are functions of classical 
phase space variables and operator product is replaced with star product [18]. The WWGM 
method works well for observables possessing a classical limit. However, it is not clear how 
it should be generalized to embrace spin degrees of freedom. Spin may be incorporated into 
classical mechanics considering the semiclassical approximation as well as the nonrelativistic 
limit of the Dirac hamiltonian. The latter is given in terms of operator valued matrices. Hence, 
we consider observables which are matrices whose elements are functions of classical phase space 
variables but depend on h. Dynamical equations of these matrix-valued symbols will be given 
by a semiclassical bracket acquired from the Moyal bracket. 

When there are some different types of gauge fields, they can be incorporated into the 
hamiltonian formalism by considering an enlarged lagrangian system which leads to second 
class constraints. Indeed, to embed Berry gauge fields in the semiclassical scheme a constrained 
hamiltonian dynamics will be presented starting from an appropriate matrix valued lagrangian. 
We adopt the Dirac formulation of constrained hamiltonian dynamics replacing the Poisson 
bracket with the proposed semiclassical bracket. This furnishes us with a systematic formulation 
of dynamics when non-abelian Berry gauge fields are present. 

Once this formulation of matrix valued observables coupled to gauge fields is accomplished 
we can employ it to investigate dynamical properties of diverse spin systems. The semiclassical 
dynamics of Bloch electrons in the adiabatic approximation where interband interactions are 
neglected was discussed in [8],[II]. We study the same problem within our approach. We 
achieved the correct phase space measure and noncommutativity of phase space variables. 
Unitary transformations which generate Berry gauge fields are also considered. We derived the 
equations of motion which can be used in topological spin transport [12]. In [15] was shown 
that when an electric field is applied to an electron, a transverse force on the spin current 
occurs. This force resulted in the Heisenberg equation of motion of velocity, considering the 
nonrelativistic limit of the Dirac hamiltonian. We show that it can easily be derived within 
our formulation. On the other hand the nonrelativistic limit of the Dirac hamiltonian can 
also be obtained by a momentum dependent Berry gauge field which generates the Thomas 
precession [20]. We calculate the force acting on an electron in the electric and magnetic fields 
in the presence of this gauge field. The same transverse force on the spin current occurs which 
depends only on the electric field. However, the terms depending both on the magnetic and 
electric fields do not concur. Experiments may settle this disagreement. The intrinsic spin 
Hall effect was envisaged in [21] analyzing the spin current due to the Rashba hamiltonian [22]. 
Investigating the Rashba Hamiltonian is still attractive, although when the vertex corrections 
are taken into account the originally proposed universal behavior of spin Hall conductivity does 
not survive [2B] (for a review see t24j). We study the Rashba spin-orbit coupling within our 
semiclassical approach to attain a very simple formulation of the spin Hall conductivity. It 
is inspired by the derivation of the Hall conductivity by demanding that the force acting on 
electrons vanishes. 

In Section 2 we present the semiclassical bracket of matrix valued observables and its basic 
properties. The semiclassical constrained hamiltonian formulation which leads to a systematic 
approach of analyzing dynamical systems with different sorts of non-abelian gauge fields is given 
in Section 3. An application of the formalism to various dynamical spin systems is considered 
in Section 4. We clarified some persisting discussions and also gave a simplistic formulation of 
the spin Hall conductivity. The obtained results and other possible applications are discussed 
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in the concluding section. 



2 Semiclassical symbols and the Moyal bracket 

Let us deal with the classical canonical variables (p^, a;^) corresponding to the quantum phase 
space (p^, x^); fi = 1, ■ ■ ■ , M. In the WWGM method of quantization one considers the symbol 
map [IS] 

S(fip,x))=fip,x), (1) 



where f{p,x) is the c-number function corresponding to the operator f{p,x). 
Let the operator product of the quantum observables / and g be 

fip,x)g{p, x) = h{p, x). 

Symbol map should respect the operator product, so that we should introduce a (star) product 
satisfying 

S {j{p, x)g{p, x)^ = S (h{p, x)) = S {j{p, x)) -k S [g{p, x)) . (2) 

Obviously the symbol map as well as the *-product depend on the operator ordering adopted. 
We deal with the Weyl ordering where the associative star product is 



exp 



ih I d d d d 



2 \ dxi^ dp^ dpf^ dx^ 



(3) 



The arrows on the derivatives indicate the direction in which they should be applied. We 
adopt the Einstein convention, hence the repeated indices are summed over. To imitate the 
commutator of operators, we define the Moyal bracket of two arbitrary observables f{p, x) and 
g{p,x) as 

[f{p,x),g{p,x)]^ = f{p,x)-kg{p,x) - g{p,x)-k f{p,x). (4) 
Hence, the classical phase space variables satisfy the Moyal bracket 

[p„x1. = -zm;:, (5) 

analogous to the canonical commutation relations. The classical limit of the Moyal bracket (jlj) 
is the Poisson bracket: 

— i df Og Of Og 

lim-[/(p,x),,(p,x)]. = {fip,x),gip,x)} ^ —— - — — . (6) 

When one considers the Dirac hamiltonian or higher spin formalisms, it is still possible 
to define a symbol map. Now observables are matrices which take values in classical phase 
space[25].[26j. The Moyal bracket of the matrices Mabip^x) and Nab{p,x) can be defined as 

{[M{p, x), N{p, x)l)^, = Mac{p, x) k N,bip, x) - Nacip, x) ^ M,b(p, x). (7) 

However, the classical limit (E]) of ([7]), in addition to the Poisson brackets of matrices, yields 
a commutator of matrices which is singular. Generally, the observables in a block-diagonal 
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form are taken into acoount for getting rid of the matrix commutator. As far as observables 
possessing a direct classical interpretation are considered, this restriction seems necessary for a 
semi classical study [26]. Indeed, we will relax this condition. When interactions are considered, 
the nonrelativistic limit of the Dirac hamiltonian may include the spin. Then there will be terms 
depending on h whose classical limit is not direct. We would like to study the semiclassical spin 
dynamics. Thus, although we deal with the classical phase space we let the symbols depend on 
h. Therefore, instead of the classical limit ([6]) we deal with the limit obtained from the Moyal 
bracket ([7]) by retaining the terms up to h: 

{M{p,x),N{p,x)}c = ^[M,iV] + ^{M{p,x),N{p,x)} - ^{N{p, x), M{p, x)}. (8) 

We would like to emphasize that the first term is the commutator of matrices, it is not the 
quantum mechanical one. Hence, it is not an attempt to combine the quantum commutator 
and the Poisson bracketH. Although we keep terms up to h order in the Moyal bracket ([7]), 
remember that M and N can depend on h. In fact, ([8]) is an expansion in powers of h where 
only the first two lowest nonvanishing terms are retained. 

Multiplication of observables is still given by the star product Q. Hence the Jacobi identity 
which should be satisfied is given by 

{M, {N, L}c}. + {N, {L, M}c}. + {L, {M, N}c}. = 
[M, {N, L}] - [M, {L, N}] + {M, [N, L]} - {[N, L], M} - i[M, [N, L]] 
+ (cyclic permutations of M, A^, L) + 0{h) = 0. 

In fact one can show that it is fulfilled up to h order. Moreover, one can observe that the 
Leibniz rule defined as 

{Mi.N,L}c = {M,L}c^N + M*{N,L}c (9) 

is also satisfied at the h order. 

To define semiclassical dynamical equations we propose to replace the Poisson bracket in 
classical dynamical equations with the semiclassical bracket ([H]). Let the symbol of the Dirac 
hamiltonian or its nonrelativistic approximation be the matrix H{p,x). Thus we consistently 
establish 

Mip,x) = {Mip,x),Hip,x)}c, (10) 

as the time evolution of the semiclassical observable M{p, x). It is worth to recalling that, as it 
is elucidated above, in this equation of motion one retains the lowest two nonvanishing terms 
in h. 

3 A semiclassical constrained hamiltonian system 

When a classical system is described with a Lagrangian, the definition of canonical momenta can 
yield some relations between coordinates and momenta which are called primary constraints. 
Preserving these constraints in time may produce some other constraints [19j. Once all the 

^For the attempts of combining the Poisson and quantum brackets, see [27] and the references given therein. 
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constraints are derived each one can be classified as first or second class due to their Poisson 
bracket relations. A method of treating second class constraints is to introduce Dirac brackets 
which effectively set the constraints equal to zero. We will consider a constrained hamiltonian 
system utilizing the semi classical bracket ([H]) and the dynamical equation (ITUi) . 
Let us consider the first order lagrangian which is a x matrix: 

C = r'' (^hy^ + pA^{r,y)+r]a^{r,y)^ - (^hr^ - ^B^{r,y)^ -no{r,y). (11) 

Here a = l,---,n, and for the nonrelativistic case the dot over the variables indicates the 
derivative with respect to time t and for the relativistic formalism it is the derivative with 
respect to an evolution parameter r. p, ^, and rj are coupling constants corresponding to the 
gauge fields A, B, and a which are N x N matrices. / is the unit matrix. Observe that in (ITTIl 
generally one cannot get rid of A, B, and a terms by redefining the coordinates or y^. The 
definition of canonical momenta 

n" = ^, = ^ 

leads to vanishing of the relations 

i)^"" = (n^ - ^?/°)J - p^" - r/a", (12) 

^|J'" = (^° + V)/-es^ (13) 

which are called primary constraints. In terms of the canonical hamiltonian Tio we need to 
introduce the extended hamiltonian 

He = no + A:C, (14) 

where A" are Lagrange multipliers and z = 1, 2. To employ the semiclassical approach of Section 
2, we identify the canonical variables as = (n°,n^) and = {ya,ra)- The semiclassical 
brackets between the constraints can be shown to be 



{^a,V'/3}c = -9ap + i^-p^-Tl^-^[A^,Bp\-—[a^,Bpl 



where 

gap is the flat metric and field strengths are defined as 



dap da^ irj 

dAp dAa ipr. , . 

= -g^-^-j[^a,Ap], (16) 

_ dBp dB^ li 



5 



Therefore, constraints (fT2l) . (fT3l) are second class and the condition of preserving them in time 

{V^^,He}c~0, (18) 

where ~ indicates that the equahty is vahd up to vanishing of constraints, will determine A^. 
In fact, in terms of 

C% = {C, <}c; C^^^'CZ'f = 5^Jl,, (19) 
one can show that (fT8l) leads to 

K = -{V;.Ho]cC-^r^. (20) 

To set effectively the second class constraints f|T2|) . f|T3l) . equal to zero, we introduce the 
semiclassical Dirac bracket 

{M,N}cD = {M,N}c - {M,r}cC;Xr\N}c. (21) 

Now, in dynamical equations the semiclassical bracket of observables ([8]) should be substituted 
with the semiclassical Dirac bracket (!2Tll . Observe that the coordinates satisfy 



{r-,r^}^^ = Cr/"^ (22) 
{y'',y^}cD = 0^2"^, (23) 
{r",/}ci. = 0^2'"'' ■ (24) 

We omitted the unit matrix I on the left hand sides. Obviously, Cj^^"'^ = —C2i°'^ = gap + ■ ■ ■ , 
thus one should consider as coordinates and ya as the corresponding momenta. 

The equation of motion of an observable 0{r, y) is given with the extended hamiltonian as 

dir,y) = {Oir,y),H,}c, (25) 
in accord with the constrained dynamical systems. Plugging the solution fl2U]) into f[T^ yields 

The inverse matrix elements C^l"^ will be obtained as a power series in the coupling constants 
p, ^ which may be identified with h. Then in the equation of motion (125!) we will retain the 
lowest two nonvanishing terms in h. 

4 Spin dynamics 

Within the formulation of the previous section we will focus on some different approaches of 
studying semiclassical dynamics of electrons in terms of Berry gauge fields. Before considering 
specific systems let us present the general formulation where Oq, = aa{r) is an abelian gauge 
field and the coupling constants are rj = e/c, ^ = h, p = —h. In our notation e < for an 
electron. The matrix C^^ defined in (fT9|) reads 
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where 

^"^ = f^ + ^ + ^[^"'^/^]- (27) 

Obviously Map does not possess any symmetry or antisymmetry with respect to the indices, so 
that one should distinguish Map from M^a- The inverse of (l26l) can be calculated at the first 
order in h as 

^iLp = ^Gal3, (28) 

Cuap = 9aP + hMpa--h{Gf)ap, (29) 
= -9ap-hMap + -^h{fG)af3, (30) 

C22aP = -Jap-hFap + ^iMf)ap-^iMf)pa-^ifGf)ap. (31) 

The equations of motion of the phase space variables can be obtained as 



^ + ^t-^^' ^{W~ 'f^'' (^"' ^ ~ V^G/)"^) , (32) 

- ^[S^, Ho]^ (^r^ - nF"'^ + ^(M/)"^ - ^(M/)^" - ^(/G/)"^)(,33) 



dHo 



+ 



at the first order in fi, employing definition ([25 



4.1 Phase space measure 

In [H],[Il] the Berry phase emerges because of keeping only lower band effects in studying 
the semi classical dynamics of Bloch electrons. To understand this formalism let = aa{r) 
be the electromagnetic gauge field with the coupling constant rj = e/c and the Berry gauge 
fields be Aa = 0, and Ba = Ba{y) with ^ = h. Although, in [8], [11] only the abelian gauge 
field was considered we let Ba be non-abelian. Hence, the matrix C is given as in (!26|) with 
= Map = 0, and at the first order in h the following semiclassical Dirac brackets result 

{ra,rp}cD 
{ra,yp}cD 

{ya,rp}cD 

{ya,yp}cD 

Similar relations were obtained in studying the electromagnetic interactions of anyons. 
Now, the equations of motion of Tq, and ya can be straightforwardly derived from ( l32l) and ( l33l) . 
respectively. 



hGap, 

9aP —{Gf)ap, 

—QaP + —ifG)aP, 

6 C^fl 

-faP :rifGf)ap- 

C 



(34) 
(35) 

(36) 

(37) 
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Adopting the formalism of the usual constrained hamiltonian systems [29], [30], the semi- 
classical phase space volume element in the presence of second class constraints is given by 

YldW^dW^dyadr^ J det ^^^C5{tpi)5{^p2)- 

\ a / 

After eliminating 11,. and IIj^ by employing constraints (|T2l) . (|T3l) and using ( |26l) with = 
Ma/B = 0, the phase space volume element becomes 



Hdy.drJ det'/'C = (j[dy.drA (^1 - ll^^ . 



(38) 



This is the phase space volume element discussed in [8],[TT]. Although in a different context 
in [S] the role of second class constraints in defining the phase space volume element f l55]) was 
noted. 



4.2 Unitary transformations 

The nonrelativistic approximation of the Dirac hamiltonian interacting with external fields can 
be obtained in terms of the Foldy-Wouthuysen unitary transformation U. In [T] and [TH] Foldy- 
Wouthuysen transformations were engaged to introduce Berry gauge fields. In [13] a projector 
on the positive energy space V is employed to define 

= VU^. (39) 

Here yi are the components of the three vector y and the flat metric is gij = Sij; i,j = 1, 2, 3. 
The Berry gauge field can be shown to bc[13j 

' 2 + mc'E,) ' ^ ^ 

where = {y ■ y)c^ + m^c^ and (Xj are the Pauli matrices. This is a non-abelian gauge field. 
Using (go]) in the general approach (l22])-([2l]),([28D-([3lD with ^ = 0, a = and ^ = /i, yields 

{yuyj]cD = 0, (41) 

{yhTjjcD = -Sij, (42) 

{r.,r,}cn = [rna + ) ■ (43) 



These coincide with the noncommutativity relations obtained in [13] . 

On the other hand in [12] a unitary transformation U = U{r, y) which diagonalizes the initial 
matrix valued hamiltonian was introduced. Generally U = U{r, y) depends on all phase space 
variables. Hence one can define the gauge fields which are non-abelian as Af = —U^^, Bf = 
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f/^^, with ^ = h, p = —h. Because of being pure gauge fields their field strengths vanish: 
= 0, Gfj = 0. However, in the adiabatic approximation one deals with 



- diag [U:^ j , Br - diag [U:^ j . (44) 
Though these are abelian gauge fields, their field strengths 

do no longer vanish. 

The equations of motion of the phase space variables can be read directly from fl32l) and 
(1331) as 



' dvj dyj 
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where /jj is the electromagnetic field strength, rj = e/c and 

n /I {ad) 

In terms of these equations of motion one can study topological spin transport. 
4.3 Transverse spin force 

We would like to discuss the spin dependent dynamics obtained in [15] within our approach 
by using the equations of motion (15^ and ([2SD- For this purpose we choose the canonical 
hamiltonian to be 

Ho = ^y^ + V + pB(T-B, (45) 
2m 

where Hb = ~eh/2mc and at the first order in h we take Ves = V(r) + ^^^2 ^ ^ 

Bi = ^eijkf-^^ is the external magnetic field and ai are the Pauli matrices. In accord with [15] 
we let i3 = and the other Berry connection be 

4mc^ dr^ ^ ^ 

Moreover, we set rj = e/c, and p = —h. We would like to emphasize that and yi are the 
coordinates and momenta in the restricted phase space. When we plug the gauge field (H6|) and 
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the canonical hamiltonian (H5!) into the equations of motion (!32l) . (!33|) we obtain 



Ti = -K— = — , (47) 
oyi m 

The force can directly be read from (H8i) in terms of the velocity v = r given by ( H71) as 

d e h f ■ d'^V ■ d'^V \ 

Ti = yi= mri = -— {V + ^ibo- ■ B) + -eijkVjBk + ^ eijkcr^vi—— h ejkia^v^ 

OTi c Amc^ \ oriOTk oriOTi J 

I^B ( Q ^ dV_ dV_ 

2mc2 \ ' ' dri dri J Sm'^c^ ^tjkO'i '^j ■ 

Indeed, this is the force obtained in [15j. The last term is the transverse spin force on the spin 
current quadratic in the electric field. 

Equations of motion following from the nonrelativistic approximation of the Dirac hamilto- 
nian can be derived in electrodynamics employing the Thomas precession [3Tj without referring 
to the Dirac hamiltonian. The relation between the nonrelativistic limit and the Thomas pre- 
cession was clarified in |20] by showing that the latter should be considered as a Berry phase 
when the external electric potential is smooth. The related gauge field can be obtained in the 
nonrelativistic limit from (HOl) . Hence, to obtain the force acting on an electron in the external 
electric and magnetic fields it should be possible to consider either the gauge field A given 
in (H6l) or the gauge field B obtained from fHOj) in the nonrelativistic limit: Let ^ = and 
deal with the electrodynamic gauge field ai(r), rj = e/c and the nonrelativistic limit of the 
non-abelian gauge field (HOl) 

1 

Field strength of this gauge field can be calculated to be 



^« = Tz::^^^lky'^'■ (50) 



-1 k 1 



where we used ^ = h. 

By ignoring the fG terms in (l32i) and (l33i) the equations of motion are 

= + —e^mBk, (52) 

r, = yi + hG,,^ + -^(B,{y.<T)-a,iB.y)). (53) 

Now, by keeping the terms linear in the velocity Vi one can show that 



mri = m{ri, He}c = + l^^o- ■ B) + ^— ^e^jfeCi^; ^^^^^^ + -eijkVjBk 



fiB f ^dV ^ dV\ h dV dV 

2mc^ \^ ' dri dri J Sni^d^ ^tjk^i '^j ■ 
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Up to some d'^V / dridrj terms this coincides with (H9|) . In fact, the latter approach is valid for 
potentials changing slowly. However, neglecting the fG terms in (!32l) . (l33|l is not justified, due 
to the fact that they may give contributions of the /is/mc^ order to the force. Indeed, retaining 
the fG terms in and using /i^ = —eh/2mc, the equations of motion of the h order are 

^iV + ^^BCT■B) ^ e f^B f dV dV\ 

^iJkVjBk - -3 -TTj^jB' - ^iBjjrr . (55) 



£ + ^^-S" - i^'^y ■ + ^^^^ ■ - ^y^^^ ■ "0 • ^^^^ 



Hence, the force linear in velocity becomes 

mVi = (V/ + yUsCr ■ B) + :j(^ijk(^-'vi^-^ h - 

ovi Imc^ oriOVk c 



+ J^(saB — + Ba — -4—(Ba)]+^—e a—v— (57) 
2mc^ \ ^ dri dri dvi ) Snfid^^^^^^^ dri^^ dr^' 



The last term which is the transverse spin force on the spin current results to be the samqj. 
However, the terms which depend on both the electric and magnetic fields are in dispute with 
(1491) . This discrepancy between the two nonrelativistic approximation schemes can be settled 
by experiments. 

4.4 The Spin Hall Effect 

Electrons constrained to move on a plane in the presence of a uniform external magnetic field 
perpendicular to the plane deviate and produce an electric field which is perpendicular to both 
the initial direction of the current and the magnetic field. This is the Hall effect which manifests 
itself as the Hall conductivity. We would like to present a derivation of the Hall conductivity 
which will inspire a simple formulation of the intrinsic spin Hall effect utilizing our semiclassical 
approach. To this aim let us deal with the hamiltonian 

H, = ^y^ + V{n,r2), (58) 
where the scalar potential is given in terms of the uniform electric field components Ei as 

V{ri,r2) = -eEiTi - ei?2'r2- (59) 

In order to constrain the electron to move on rir2-plane we set = 0. We consider the vanishing 
Berry gauge fields A = 0,B = and let there be a uniform magnetic field in direction: 

/i2 = B. (60) 

The related coupling constant is rj = e/c. The equations of motion following from ( l32l) and ( l33l) 
are 

r^ = ^, (61) 
m 

yi = eE^ + —y2, (62) 
mc 

1/2 = eE2 - —1/1, (63) 
mc 



^In [32] it was claimed that this method leads to a transverse force in conflict with [TSl. 
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The force acting on electron can be read from (l6T]) - (!63!) . in terms of the velocity v = r, as 



eB 

JFi = mfi = yi = eEi + —^2, (64) 

J^2 = = m = e-E'2 - —Vi. (65) 

c 

Till now we have considered single particle dynamics. To connect it to a system of electrons 
let us introduce the density of electrons k. Thus, the electric current is defined by 

j = env. (66) 

We demand that the net force acting on electrons vanish Ti = 0, so that the electrons move 
without deflection (see e.g. [33]). We can solve this condition for the velocity and plug it into 
( !66l) . which yields the electric current 



J2 J [(TH J [ E2 



(67) 



where 

ecK 

^^ = -5- 

is the Hall conductivity. 

The intrinsic spin Hall effect is envisaged in |21] in terms of the Rashba spin-orbit coupling [22] 
By generalizing the Hall effect formulation we can introduce a simple method of acquiring the 
spin Hall effect conductivity employing the Rashba spin-orbit coupling. The hamiltonian is 
still given by flSHjl with y-^ = 0. However, there is no magnetic field: = 0. To consider the 
linear Rashba spin-orbit coupling we set B = and define 

Ai = eijkUjel. (68) 

Here is the unit vector in the third direction el = 6k3 and o", are the Pauli matrices. Moreover, 
in the original formulation (fTTj) we should take p = —am/h, where a is the Rashba coupling 
constant |22j. 

The related field strength can be calculated as 

Fij = -j[Ai, Aj] = y (73e,,fce^. (69) 

The equations of motion of the canonical variables are 

= ^, (70) 
m 

m = + (71) 

Hence, the force acting on the particle is 

2p2 

JF. = mri = eEi + —cr^eijkelvj. (72) 
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Imitating the formulation of the Hall effect we set jFj = 0, in order to have a motion without 
deflection. This condition is solved for the velocity as 



,^-^E2,vi = -0E2, (73) 

vl = -j^E^,v^ = j^E,. (74) 

The arrows t and | correspond, respectively, to the positive and negative eigenvalue of (73. It 
is natural to define the spin current as 

= (^n^v^ — n^v^) , (75) 

where and denote the concentrations of states with spins along the and — directions. 
Employing ([73]), ([74]) in ([75]) yields 



where 



r = asHe' X E, (76) 

is the spin Hall conductivity. In this simplistic approach the total concentration of states 
n = (n^ + n^^ is an input which should be given by other means. Although its calculation is 
beyond the scope of this work, for having an insight let 

n = snl^ (78) 

where s is a constant and ng^) is the concentration of states occupying the lower energy state 
of the Rashba hamiltonianf 



2 2 
am 



n. 



2^ " vrn^ • 

Using (|78ll in (|77l) leads to the spin Hall conductivity 

CTSH = (79) 

This agrees with the universal behavior obtained in pi] for s = 1/2. However, when the vertex 
corrections are taken into account it is known that this universal behavior does not survive[23j, 
as far as the linear Rashba coupling is considered. The vertex corrections were calculated 
employing Green functions within the Born approximation. Hence, it is not clear how one 
can incorporate the vertex corrections into our semiclassical scheme. To cure the defects of 
the linear theory, it would be useful to study the Rashba couplings which are higher orders 
in momenta (see |21j and the references therein). Although we will not discuss it here, our 
semiclassical approach can be used to investigate higher order generalizations of the Rashba 
spin-orbit coupling. 
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5 Discussions 



Semiclassical limit designated as the bracket ([8]) can be utilized to study diverse dynamical 
problems where spin degrees of freedom are not ignored. Hence, instead of dealing with wave 
packets one can consider single particle interpretation of semiclassical dynamics of spin depen- 
dent systems. 

It can be shown that the constrained hamiltonian system which we presented here is suitable 
to investigate properties of some topological quantum phases. Moreover, as we will present in 
a future work it constitutes a new gauge invariant method of studying dynamical systems in 
noncommutative spaces. 

Any model concerning spin dynamics utilizing Berry gauge fields which give rise to non- 
commutativity of coordinates and/or momenta can be studied in terms of the semiclassical 
approach presented here. We focused on some recent formalisms where some of persisting dis- 
cussions can be clarified. The results which we derived are valid up to the first order in h. 
However, in this formulation keeping the track of the higher orders is possible. When higher 
order h corrections are considered there may be some different sources: Gauge fields which we 
consider may depend on higher h, the limit of the Moyal bracket will have another term and 
inversion of the matrix C^^ may lead to some higher h terms. 

Obviously, the formalism of the spin hall conductivity which we reported here should be 
elaborated. Nevertheless, due to its resemblance with the Hall effect and simplicity, it may be 
profitable to predict some basic properties of the spin Hall effect. 
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